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Abstract 

A less studied numerical characteristic of periodic orbits of area preserving twist 
maps of the annulus is the twist or torsion number, called initially the amount 
of rotation f23]. It measures the average rotation of tangent vectors under the 
action of the derivative of the map along that orbit, and characterizes the degree of 
complexity of the dynamics. 

The aim of this paper is to give new insights into the definition and properties of 
the twist number, and to relate its range to the order properties of periodic orbits. 
We derive an algorithm to deduce the exact value or a demi-unit interval containing 
the exact value of the twist number. 

We prove that at a period-doubling bifurcation threshold of a mini-maximizing 
periodic orbit, the new born doubly periodic orbit has the absolute twist number 
larger than the absolute twist of the original orbit after bifurcation. We give exam- 
ples of periodic orbits having large absolute twist number, that are badly ordered, 
and illustrate how characterization of these orbits only by their residue can lead to 
incorrect results. 

In connection to the study of the twist number of periodic orbits of standard-like 
maps we introduce a new tool, called 1-cone function. We prove that the location of 
minima of this function with respect to the vertical symmetry lines of a standard- 
like map encodes a valuable information on the symmetric periodic orbits and their 
twist number. 



1 Introduction 

The study of the dynamics of area preserving positive twist maps of the annulus is mainly 
concerned with the characterization of its invariant sets, and the dynamical behaviour 
of a map restricted to such sets. Among invariant sets, rotational periodic orbits have 
been thoroughly studied and classified according to their linear stability, extremal type, 
order properties (see [21] for a survey). A numerical characteristic associated to a periodic 
orbit is the rotation number, which measures the average rotation of the orbit around the 
annulus. In [23j Mather defined also the amount of rotation, which is called twist number 
in p] or torsion number in [lO]. The twist number of a periodic orbit characterizes 
the average rotation of tangent vectors under the action of the tangent map along the 
map's orbit. Mather related the twist number with the Morse index of the corresponding 
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critical sequence (configuration), and Angenent [2] proved that in the space of {p,q)- 
sequences a critical point of the W^pg-action, corresponding to a periodic orbit of twist 
number greater than is connected by the negative gradient flow of the action, through a 
heterochnic connection, with a sequence corresponding to an orbit of zero twist number. 
Using topological arguments, Crovisier [10] proved the existence of orbits of zero-torsion 
(twist) number for any real number in the rotation number set of a twist map of the 
annulus (not necessarily area preserving). 

In order to detect new classes of dynamical systems exhibiting periodic orbits of non- 
zero twist (torsion), Beguin and Boubaker gave |5] conditions ensuring that some area 
preserving diffeomorphisms of the disc D^, and particular diffeomorphisms of the torus 

exhibit such orbits. 

In Mather's and Angenent 's definitions of the twist number this a positive quantity. 
However its natural definition leads to a signed number [ID]. The twist number of periodic 
orbits of positive twist maps is zero or negative, while for those of negative twist maps it 
is zero or positive. In our approach we keep its natural sign and call absolute twist, the 
absolute value of the twist number. 

During almost 30 years since the definition of the twist number was given, no periodic 
orbit of absolute twist number greater than 1/2 was detected in the dynamics of classical 
twist maps (standard map, Fermi-Ulam map, etc). That is why we wonder under what 
conditions can such orbits appear. 

From Aubry-Mather theory [3], [22] it is known that non-degenerate minimizing and 
associated mini-maximizing [p, 5')-type periodic orbits are well ordered and the minimiz- 
ing orbit has zero twist number. Thus it is natural to investigate whether (p, g)-periodic 
orbits that are not given by this theory can be well-ordered or not and how order prop- 
erties are related to the value of their twist number. 

Our starting point was not however the study of the twist number. In an attempt to 
characterize the dynamical behaviour of twist maps after the breakdown of the last KAM 
invariant torus (the case of standard map) or in a half annulus where no invariant circle 
exist (the case of Fermi-Ulam map [TS] , or tokamap ^ ) , we identified in the phase space 
of such maps, special regions where no minimizing periodic orbits can land. Instead we 
noticed that any periodic orbit having at least two points in such a region is badly ordered 
and typically has absolute twist greater than 1/2. 

Thus we were led to a deeper study of the twist number of periodic orbits, and their 
properties. In this paper we complement the results on twist number reported in [23] . 
[2], and show that periodic orbits of non-zero twist numbers are typically born through a 
period doubling bifurcation. We give examples of periodic orbits that have large absolute 
twist number. As far as we know no such orbits were identified before in the study of 
the most known twist maps, standard-like maps. Some of these orbits are ordered and 
other are unordered. In order to explain why some twist maps can exhibit a sequence of 
bifurcations of a periodic orbit leading to increasing the absolute twist number up to its 
maximal value we introduce two subclasses of standard-like maps, USF maps and TSF 
maps. 

The paper is organized as follows. In Section |5] we give the basic properties of twist 
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area preserving diffeomorphisms of the annulus, relevant to our study. In order to derive 
in Section H] sufficient conditions that favor the existence of ordered periodic orbits of 
large absolute twist number, we set up in Section |3] a framework of study, defining the 
class of twist maps exhibiting the so called strong folding property, and a function, called 
1-cone function. 

The 1-cone function is defined on the phase-space of a twist map, and takes negative 
values within the region where the map exhibits strong folding property. We prove that 
the restriction of this function to a periodic orbit gives information on the eigenvalues of 
the Hessian matrix associated to that orbit (Lemma 13. ip . Analyzing the 1-cone function 
associated to a standard-like map, we show that such maps can have either a connected 
strong folding region including one of the vertical symmetry lines or a two-component 
strong folding region including both vertical symmetry lines. These results will be ex- 
ploited in Section 5 to explain the variation of the twist number. 

In Section m we revisit the definition and properties of the twist number of a periodic 
orbit based on the structure of the universal covering group of the group SL{2,'R) (a 
system of coordinates on this group allowing to decipher its topology is presented in [B|. 
The twist number is defined as the translation number of a circle map induced by the 
monodromy matrix associated to the periodic orbit. We review in|X]the properties of the 
translation number of an orientation preserving homeomorphism of the unit circle and 
point out the particularities of the translation number of circle homeomorphisms induced 
by matrices in SL{2,'R). 

Theorem 4.2 gives the relationship between the twist number value of a {p, g)-periodic 
orbit, and the position of the real number with respect to the sequence of interlaced 
eigenvalues of the Hessian matrix Hg, associated to the corresponding {p, g)-sequence, 
and of a symmetric matrix derived from Hg. This theorem complements results from [23] 
and [2]. 

Based on this theorem we derive an algorithm to deduce the exact value of the twist 
number or a demi-unit interval that contains the twist number. 

The main result in Section H] is the Proposition 14. H which shows that periodic orbits 
of large absolute twist number are born through a period doubling bifurcation. More 
precisely it states that if at some threshold, a period doubling bifurcation of a mini- 
maximizing {p, g)-periodic orbit occurs, with transition from elliptic to inverse hyperbolic 
orbit, then the new born 2g-periodic orbit is elliptic, having the twist number within the 
interval (-1, -1/2). 

In Section |5] we give examples of periodic orbits that have large absolute twist number 
and are unordered. A natural question is whether a positive twist map can also exhibit 
ordered (p, g)-periodic orbits {p,q relative prime integers), of twist number less than 
—1/2 (we note that periodic orbits of large absolute twist number, born through a period 
doubling bifurcation are of type (2p, 2g)). We give a positive answer to this question, 
giving an example of three-harmonic standard map having an ordered (1, 2)-periodic 
orbit, of twist number r = — 1. 

In order to illustrate the different range of the twist number of a periodic orbit in the 
three-harmonic standard map in comparison to that of a periodic orbit of the same type 
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of the standard map, we deduce (Proposition 15.11) the bifurcations that a periodic orbit 
of a family of standard-hke maps, with a uni-component strong folding region, undergo 
necessarily. 

The concrete examples of periodic orbits of large absolute twist number, given in this 
section, illustrate that the classical method used for more than 30 years in classification 
of periodic orbits can lead to an erroneous conclusion. 



2 Background on twist maps 

We recall basic properties of twist maps, relevant to our approach. For a detailed presen- 
tation of classical results concerning dynamics of this type of maps, the reader is referred 

to na, m- 

Let S"*^ = M/Z be the unit circle, A = S-"^ x M, the infinite annulus , and tt : — )■ A, 
the covering projection, 7i{x,y) = (x mod l,y). 

We consider a C"^-diffeomorphism F : — i- M?, F{x,y) = {x',y'), satisfying the 
following properties: 

i) F is exact area preserving map, isotopic to the identity; 

ii) F{x + l,y) = y) + (1, 0), for all (x, y) G M^; 

dFi 

iii) F has uniform positive twist property, i.e. — — > c > 0. 

dy 

The map F defines a C^-diffeomorphism, / : A — )■ A, such that vr o F = / o tt. Both 
maps / and F are called area preserving positive twist maps or simply twist maps. F is 
a lift of /. 

In the sequel we will switch from / to F, and conversely, without comment. 
The motion of a point around the annulus is characterized by its rotation number. 
The orbit of a point z e A has a rotation number if there exists the limit: 

p = lim — — -, (1) 

71— >oo n 

where {x,y) G is a lift of z, and {xn,yn) = F"'{x,y). The rotation number does not 
depend on the chosen point z on the orbit or the lift {x,y). For different lifts of the map 
/, the corresponding rotation numbers differ by an integer. 

Let p, q be two relative prime integers, g > 0. A {p, g)-type orbit of the twist map, F, 
is an orbit {xn,yn) = -^""(3^0)1/0); n E Z, such that Xn+q = Xn + p, for every n G Z. The 
vr-projection of such an orbit onto the annulus is called {p, q)- periodic orbit or simply 
g-periodic orbit of the map /, and its rotation number is p/q. 

The linear stability of a {p, g)-periodic orbit, 0{z), is characterized by the multipliers 
of the derivative, D^f^ G SL[2,M.), which in turn are determined by the trace, ii{Dzf'^)- 
Instead of trace one can also use the residue [H] to classify the g-periodic orbits. Namely, 
the residue is defined by: 

^ ^ 2-trace(Z}jg) 
where q is the least period of the orbit. 
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The dynamics of a twist map, F, has a variational formulation, that is the F-orbits 
are associated to critical points of an action, and vice-versa. We discuss only the case of 
{p, g)-type orbits. 

More precisely, F admits a C^-generating function /i : — M (unique up to additive 
constants), such that hi2{x,x') = dh/dxdx' < on M^, and F{x,y) = {x',y') iS y = 
—hi{x, x'), y' = h2{x, x') {hi = dh/dx, h2 = dh/dx'). 

The generating function defines an action Wpq, on the space Xpq, of (p, g)-type se- 
quences of real numbers, x = (x„)„g2, i.e. sequences satisfying the property, Xn+q = 
Xn + p, for all n G Z. The space Xpq is identified with the affine subspace of M'^"'"^ = 
{(xo, xi, . . . , Xq)}, of equation Xq = xq + p, and the action is defined by: 

9-1 

Wpq{x) = '^h{Xk,Xk+l) 
k=0 

There is a one to one correspondence between critical points, x = (x^,), of the Wpq-action, 
and (p, q')-type orbits, (x„, of the positive twist map, F: 

{Xn, Vn) = {Xn, -hi{Xn, X„+l)), V 72 G N (3) 

By Aubry-Mather theory [3], [22], for each pair of relative prime integers (p, q), q > 0, 
a twist area preserving map has at least two periodic orbits of type (p, q). One corresponds 
to a non-degenerate minimizing (p, g)-sequence of the action, and by abuse of language 
it is also called in the following, minimizing {p, q')-type orbit. Associated to such an 
orbit there is a second one, corresponding to a mini-maximizing sequence, and called 
mini-maximizing {p, (3')-type periodic orbit. 

Let Hq be the Hessian matrix of the action, Wpq, at a critical point x = The 
signature (the number of negative and positive eigenvalues) of the Hessian Hq, at a non- 
degenerate minimizing sequence is (0,g), while at the corresponding mini-maximizing 
sequence it is (1, q — !)■ The number of negative eigenvalues is called the Morse index of 
the critical sequence. 

A salient property of a non-degenerate minimizing {p, g)-periodic orbit of a positive 
twist map, /, is that it is well-ordered. 

An invariant set M of a positive twist map, /, is well ordered if for every {x, y), {x', y')' G 
7r~^(M), we have x < x' \S. Fi{x,y) < Fi{x',y'), where Fi is the first component of the 
lift F. 

A well-ordered {p, g)-orbit is also called Birkhoff orbit, while a badly-ordered orbit is 
called non-Birkhoff. 

A {p, g)-periodic orbit of the twist map / is well ordered iff the coresponding (p, g)- 
sequence x is cyclically ordered, i.e. 

X < TjjX or TjjX < X, V z, j G Z, (4) 

where Tij is the translation defined by: 

(rijx)fc = Xk+i +j, V X = (xfe) G M^, i, j,k eZ 
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By Aubry-Mather theory Wpg-minimizing sequences are cyclically ordered. Moreover, 
the mini-maximizing sequences are also ordered with respect to minimizing sequences 

The space Xpg of {p, g)-sequences is partially ordered with respect to the order inher- 
ited from R^. Namely if x = (xk), y = (yk) ^ K^, then: 

x<y ^ Xk<yk,ykeZ (5) 

One also defines the following relations: 

X < y ^ X < y, but X 7^ y 

To a {p, g)-sequence x = {xk) one associates the Aubry diagram, i.e the graph of the 
piecewise affine function A : R — > M, that interpolates linearly the points {k,Xk),k G Z. 
The Aubry function associated to a minimizing {p, g)-sequence is increasing. If the Aubry 
diagram of a sequence y = (yk) do not intersect the Aubry diagram of the minimizing 
sequence x = (x^), then the two sequences are comparable, i.e. either x -< y or y -< x. 
If the two diagrams intersect transversally, then the sequence y = (yk) corresponds to an 
unordered {p, g)-periodic orbit. 

The most studied twist maps, both theoretically and numerically, are the standard- 
like maps. A standard-like map is a twist map, F^, defined by a Lagrangian generating 

function of the form h{x,x') = -{x — x'Y ~ ^y{x)-, where is a fixed 1-periodic even 
function of class C"~, r > 3, and e G M is a parameter: 

x' = X + y — eV'{x) 
y' = y-eV'{x), 

Classical standard map fl2Sl) corresponds to the potential V(x) = —- — -r cos(27ra;). 

(27r)2 

The twist map, F^, is reversible, i.e. it factorizes as = / o i?, where R and / are the 
involutions, R{x,y) = {~x,y — eV'{x)), I{x,y) = (— x + y,y). The symmetry lines, that 
are the fixed point sets Fix(i?), Fix(/), have both two components: Fq : s = 0, Fq = 0.5, 
respectively Fi : x = y/2, r\ : x = {y + l)/2. The i?-invariant orbits are called symmetric 
orbits. 

As we shall show in Section 4, the twist number of a (p, g)-periodic orbit one deduces 
analyzing the relative position with respect to of the eigenvalues of the associated 
Hessian matrix and of a companion symmetric matrix. In the following section we prove 
that diagonal entries of these matrices encode information on the position of the smallest 
eigenvalue. In order to extract such an information we introduce a new tool of study 
of a twist map, the 1-cone function, whose restriction to {p, g)-periodic orbit gives the 
diagonal entries in the two symmetric matrices. Depending on the location of the minima 
of the 1-cone function, with respect to the vertical symmetry-lines of a standard-like map 
we define two subclasses of this set of twist maps: the class of USF maps, respectively 
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TSF maps. Then in Section 5 we prove that for a family of USF maps, respectively TSF 
maps, one can predict the sequence of bifurcations that their symmetric periodic orbits 
can undergo, as well as the variation of their twist number, only from the analysis of the 
diagonal entries of the Hessian and the companion matrix. 

3 1— cone function and strong folding property of a 
twist map 

Let / be a positive or negative twist APM defined on the annulus A. The derivative of 
the map at an arbitrary point has the entries denoted as follows: 



a{z) b{z) 
c{z) d{z) 



The twist APM, /, is called (+(i)-map if d{z) > 0, \/z E A, respectively (+a)-map, if 
a{z) > 0, V^; G A. The two conditions are not equivalent. 

We study only twist maps / that are {+d)- maps, because their inverses are 
maps. 

The (+d)-condition ensures that the lines of the vertical foliation of the annulus are 
mapped to graphs of functions over the action, while the sign of a{z) illustrates how the 
map acts on the horizontal foliation of the annulus. If for a (+(i)-twist map the function 
a{z) is negative on a segment of an horizontal circle, x {y}, then the /-image of that 
circle is not a graph over x. We say that the circle is folded by the map. 

Definition 3.1 A {-\-d)-twist map is said to satisfy the folding property if there exist 
regions in the phase space where a < 0. The subset J-' = {z E A\ a{z) < 0} is called 
folding region. 

The standard map is a (+ci)-twist APM which exhibits folding property for e > 1. 
The vertical strip: 



J-e = {{x,y) I arccos(l/e) < x < — arccos(l/e)}, e > 1 

2tt 2tt 



(7) 



is its folding region. 

In the sequel we refer only to (+(i)-positivc twist maps, and all results can be mod- 
ified accordingly for (+(i)-negative twist ones. In order to reveal the properties of twist 
maps exhibiting a folding region, we use geometrical arguments, that are more intuitive 
and allow visualizations in the phase space. The same properties could be presented in 
configuration space in a more formal setting. 

We associate to each point of the phase space a pair of vectors: 



vt{z) = Df-.^,)f{e,))^{b{f-Hz)),d{f-Hz))r, 
vTi^) = Df^,^f-\-e2)^{b{z),-a{z)r, 



(8) 
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where 62 = (0,1)"^ is the vertical tangent vector at f~{z), respectively at f{z). The 
directions to which the vectors Vi{z), v^{z) are pointing illustrate how the linear maps 
Df-i(z)f, Df(^z)f~^, deviate the the upward vertical through f~^{z), respectively the down- 
ward vertical through f{z). If at a point the slopes of the two vectors, sf{z) = -{f~^{z)), 

Si{z) = —-{z), are related by > Si{z), one says that 1-cone crossing occurs at z. 

It is known that no minimizing orbit can pass through such a point |20j . 

In a folding region, since a{z) < 0, the two slopes are both positive, and 1-cone 
crossing is possible. It is easy to check that 1-cone crossing can occur in the phase space 
of a (+(i)-twist map only if the map has folding property, and the set, 

Conei = {z G A I Si{z) > s^(z)}, 

of the points where 1-cone crossing occurs is necessarily included in the folding region. 

A(+ci)-twist map with the property that the 1-cone crossing occurs in a subset of a 
folding region is said to exhibit strong folding property, and the set Conei is called strong 
folding region. 

The standard map exhibits a strong folding region for each e > 2: 

Conei(e) = {(x, y) G J-" I arccos (2/e) < x < — arccos(2/e)} (9) 

2n 2'K 

In order to show how the strong folding property influences the dynamical behaviour 
of a twist map, we define a real function 5, which associates to each point z = (x, y) the 
difference between the slopes sf{z), 

N d{x_i,y_i) a{x,y) ^ 

o{x, y) = TT r + TT {x-i,y-i) = f {x, y) 

b[x,y) 

We call this function, 1-cone function associated to the twist map. 5 takes negative values 
at the points of a strong folding region of a (-l-ci)-twist map. 

In FigJTJwe illustrate the graph of 5 over the phase space A identified to [—0.5, 0.5) x 
M (considered as a subset of the plane xOy), for different values of the perturbation 
parameter e of the standard map, while in Figl2] is illustrated the graph of the same 
function associated to the standard-like maps F^, defined by a three-harmonic potential, 

V{x) = — cos(27ra;) + ^ cos(47rx) + — cos(67ra;), (10) 
2% 4:71 Gvr 

where the parameters rji, i = 1,3, are respectively 0.18, —0.42, —0.11. 

We note that in variational setting, the function 6{x,y{x,x')) is the Hessian of the 
orbital action W{x,x') = h{x,x'), corresponding to the 1-segment of trajectory {x,y). 

Moreover, the restriction of the function 6 to the points of a {p, g)-sequence gives the 
diagonal entries of the Hessian at that sequence. Actually, if {xi,yi), i = 0,g — 1, are 
points on a (p, q')-type orbit, and 

D(x y-)F =( 
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then the Hessian expressed in (x, ?/)-coordinates is: 






(3-3 Oiq-2 Pq-2 
L ... (3q-2 &q-l 



do /3o 
/3o oil h 

■•. ■•. 







... /3, 



111 



where dj = + ajfe^ , /3j = -l/h, z = 0, g - 1. 

Lemma 3.1 Lei (xj, 6e a (p, q)-type orbit, and Hg its Hessian matrix HI]). If Xmin, ^max 
are the smallest, respectively the largest eigenvalues of Hg then the values of the 1-cone 
function, S{xi,yi), at each point of the orbit obeys the inequalities: 



Proof: The smallest and the largest eigenvalues of Hg are as follows: 



mini?(t>), 



maxR(v), 



(12) 



(13) 



where R{v) 



<i V H V > 

^ — , f 7^ 0, is the Rayleigh quotient associated to the symmetric 

< w, u > 

matrix Hg. Because Xmin < Riy) ^ Amax,Vf ^ 0, and -R(ej) = 6{xi,yi), Vi = l,g, where 
Ci is a vector of the standard basis in M'^, it follows the sequence of inequalities. 



By the Gershgorin Theorem [26] each eigenvalue, A, of a symmetric matrix A 



(ai 



belongs at least to one of the intervals 

[an - Ti {A) , an + {A)] , {A) 



E 

i=i> jV« 



"ij 1 1 



1, n 



It follows that in the case of the matrix Hg, associated to a periodic orbit, the smallest 
eigenvalue belongs to the interval: 

Xmin^[h-rk{A),5k], 6k = mh]_6{xi,yi) 

i=0,q-l 

Since the smallest value of the function 6 restricted to a periodic orbit is an upper 
bound for the smallest eigenvalue of Hg, whose sign is an indicator for the Aubry-Mather 
periodic orbits, the location of the absolute minimum of the function can give valuable 
information on extremal type of these orbits. 

Let us study the influence of location of minima of the function 6, on the dynamical 
behaviour of the standard-like maps. 

We restrict our study to standard-like maps corresponding to e > 0. Standard-like 
maps are (+(i)-twist APMs, that exhibit folding regions beyond a threshold e', which is 
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the smallest parameter e, for which the set T = {{x,y) \ a{x,y) < 0} is non-empty or 
equivalents the set {x e [0, 1) | V"{x) > 1/e} 7^ 0. 

Standard-like maps also exhibit strong folding property, beyond a parameter e* > e' 
which is the smallest parameter for which the 1-cone crossing occurs within the folding 
region: S{x,y) = 2 — eV"{x) < 0. Both folding and strong folding regions, are union of 
vertical strips, i.e. subsets of the form {{x,y) \ A < x < B,y ^ R}. 

Depending on the potential V, the corresponding standard-like map can exhibit one 
or more strong folding regions. Since V is an even periodic function, V" is also even, and 
its critical point set contains the points x = 0,x = 0.5. If at least one of these two points 
is a maximum point for V", then beyond a perturbation parameter, the corresponding 
symmetry line (i.e., Tq or/and Fq) is included within the strong folding region. Such a 
line is a line of global or local minima for the function 6{x, y), because S does not depend 
on y (see Figs. 

Remark 3.1 // one of the symmetry lines Fq, Fq, of equation x = Xq, Xq G {—0.5, 0}, is a 
line of minima for the function 6, then the restriction of the function 6 to that symmetry 
line, 6{xo,y',e), is a decreasing function with respect to e. 

We call USF map (uni-component strong folding region map) a standard-like map 
exhibiting only one connected strong folding region, and TSF (two-component strong 
folding region map), a map exhibiting a strong folding region with two connected compo- 
nents (one including the symmetry line Fq, and another, Fq). 

Standard-like maps defined by a potential V, whose second derivative has only two 
critical points is an USF map. Classical standard map fl28l) obeys this condition, as well 
as some maps having an analytic potential with infinitely many harmonics, such as the 
map: 

x' = X + y' 

, e sin27ra; , lal < 1 (14) 

y = 

27r 1 — a cos 2ttx 
For a = —0.3, for example, it is an USF map. 

The map corresponding to the three harmonic potential fIlOp is a TSF map. In Figl2] 
one can see that the function 6 associated to this map gets negative beyond a parameter, 
within two vertical strips. 

Obviously in the case of a potential given by a trigonometric polynomial the coefficients 
of the polynomial can be tuned such that the map be a USF map. 

In Section \5\ we show how the number of strong folding components influences the 
dynamical behaviour of a standard-like map, respectively the order properties of periodic 
orbits of large absolute twist number. 

In the sequel we reconsider the definition of the twist number and derive new proper- 
ties. We measure the twist number of an orbit with respect to the canonical trivialization 
of the tangent bundle of the annulus. The general framework of its study is based on the 

structure of SL{2,M.), and properties of translation number of circles maps given in the 
\M andE 
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Figure 1: The graph of the function 6 with respect to the plane xOy, associated to the 
standard map corresponding respectively to e = 0.5, 1.2, 2.3. 




Figure 2: The graph of the function 6 with respect to the plane xOy, associated to the 
standard-like map corresponding to a fixed three harmonic potential (see the text) and 
respectively to perturbation parameters e = 0.05, 0.35, 0.48 



4 The twist number revisited 

Let / be a twist area preserving (APM) of the annulus A, and zq,zi, . . . , Zq-i a g-periodic 
orbit. We identify the tangent spaces at Zi, T^. A, i = 0, g — 1, with M?. If v = is a 
tangent vector at a point z of the orbit (let us say that v G T^qA), then the corresponding 
(forward) tangent orbit is: 

{zk, Vk), where Vk+i = D^J{vk), k eN 

The amount of rotation about a periodic orbit is the average rotation of tangent vectors 
under the action of the tangent map along that orbit. In order to characterize the amount 
of rotation we interpret the derivative matrices, D^f G S'L(2,M), as circle maps (see |A|) . 
and a particular lift D^f allows then to measure how much is rotated a vector v G T^A, 
by DJ. 

To be more precise, we consider §^ as being either the unit circle, C, in or as the 
quotient group M/Z. 

Let ^ \ E SL{2, M.) be the matrix D^f, with b positive or negative depending on 

whether / has a positive twist or a negative twist. Dzf maps the vector (0, 1)"^ to (6, d)^ , 
i.e the twist map / tilts the vertical line through z G A to the right (6 > 0) or to the 
left (6 < 0). Regarding D^f as a circle map this means that the point 1/4 + Z G (the 
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north pole in = C) is displaced clockwise/anticlockwise less than 1/2. 

This property of D^f, transferred to a particular lift, amounts to set a condition on 

1 1 
its displacement map value at a; = - (x is a lift to M of - + Z G S^): 

Definition 4.1 A lift D^f of the derivative matrix, D^f , of a positive (negative) twist 
map, such that its displacement map satisfies, \E'(l/4) G (—1/2,0) (\I'(l/4) G (0,1/2)) is 
called basic lift. 

In Figj3]we illustrate the action of D^f on the unit vector v = (0, 1)^ or equivalently 
on p = 1/4 + Z G S^, when 6 > 0, as well as how the basic lift Dzf displaces the point 
X = 1/4, tt{x) = p, along M (vr : M ^ §\ 7r(x) =x + Z). 




DJix) 



Figure 3: Illustration of the action of Dzf on a unit vector v = p E S^, and of the 
corresponding basic lift, on a point x G 7r^^(p). g denotes the circle map defined by the 
matrix Dzf. 



Let Dz^^f be the basic lifts associated to matrices Dz^f G 5*17(2, M), k = 0,g — 1, 
where zq, zi, . . . , Zq^i G A are points of a {p, g)-periodic orbit, of least period q. Because 

the covering projection 11 : SL{2,M.) — )■ SL{2,M.) is a group homomorphism, each lift 

Gk = Dfk-i(^z)f o • ■ ■ o Df(^z)f o Dzf is a lift of D J^ k = l,q. We denote by the 
corresponding displacement map, and call Gk the the k-basic lift associated to the g- 
periodic orbit, k = l,q. 

Definition 4.2 The twist number of a {p,q) -periodic orbit, 0{z), is the translation num- 
ber of the associated q-basic lift, Gq: 

n— >-oo n 
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If X G M then its displacement: 



Dfi-liz)! o ■ ■ ■ o o DJXx) 



X 



colinear with 



measures the amount of rotation in a period, of a vector, v G T^A = 
u = (cos27rx, sin27rx)"^, v = fiu, /i > 0. 

The twist number defined as a translation number, is a S'L(2, R)-conjugacy invariant, 
and thus it does not depend on the point z of the orbit, to which one associates the 
g-basic lift. 

In the sequel we discuss^ the twist number of periodic orbits of positive twist maps. In 
this case the basic lift, Dzf, at an arbitrary point z G A, belongs to one of the subsets: 

T-Lq, E-i.l-i'^i C 5*17(2, M) or to the common boundary of two of them (see FigJT^. Thus 

the basic lifts Djk^^^-^, k = 0,q — 1, associated to a g-periodic orbit have the translation 
number less or equal to zero. 

The positive twist property of the map leads also to: 

Lemma 4.1 If z belongs to a q-periodic orbit, and for some 1 < k < q, the principal 
lift Gk associated to z has its displacement at 1/4, less than j/2, i.e. \E'fc(l/4) < j/2, for 
some negative integer, j , then no one of the lifts Gk+i, ■ ■ ■ ,Gq can have the displacement 
at 1/4, greater than j/2. 

In view of this Lemma and taking into account that the twist number of a g-periodic 
orbit belongs to the range of the displacement map, "^g, i.e. r G [m{Gq), M{Gq)], we can 
conclude that the twist number of a periodic orbit of a positive twist map can be zero or 
negative. 

The twist number defined as above is minus the twist number defined by Angenent [2] . 
The negative sign illustrates the clockwise sense of rotation about the periodic orbit. Its 
natural sign allows a better understanding and explanation of its properties in connection 
with the displacement map of the lift involved in its definition. 

In [23] the amount of rotation about a periodic orbit of a positive twist map is defined 
in a similar way, but in order to get a positive twist number, it was chosen for each matrix 
-D/fe(2)/, = 0, g — 1, a lift G, with the translation number t{G) G [0, 1/2). Let us argue 
that such a choice of a lift is not always possible. 

First we recall shortly the method used in [23] to establish a relationship between the 
twist number of a periodic {p, g)-type orbit, and the Morse index of the corresponding 
critical sequence of the W^pg-action. 

By a non-symplectic change of coordinates: 

(x, x') I— J- (x, y{x, x')) 

the derivative D.F'^ = Df,-^{z)F ■ ■ ■ Df{z)F ■ D^F G S'L(2,]R), at a point z 
{p, g)-orbit, is expressed as: 



M 



1 

l3q-i _ao 



1 

' P2 P2 





/3i 



1 

Ti 



(16) 
(x, y) of a 

(17) 
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The entries ai, /3i are defined by: 



(3i = hi2{xi,Xi+i), tti = hn{xi,Xi+i) + h22{xi_i,Xi), i e Z (mod g), (18) 

where h is the Lagrangian generating function of F, and hj^^, k,£ & {IjS}, are second 
order partial derivatives of h. 

Thus D^F'^ is GL+{2, M)-conjugated to the matrix M {GL+{2, R) is the group of real 
2 X 2-matrices of positive determinant), and the twist number is defined as the translation 
number of some lift of M G GL+(2, M) . 

The lifts Gb of circle maps, qb, defined by matrices B G GL"^(2,M) were considered 
in [23] as having the properties of general circle maps. The lack of information that the 
difference between the maximum and minimum displacement of a lift is less than 1/2 
led to a non concordance. 

Let us point out that the lifts Gb, B G GL+(2,M), have the same particularities as 

the lifts in 5*1/(2, M). Indeed, to each matrix B G GL+(2,M) one associates the matrix 



Thus a lift of the circle map defined by a matrix in S G GL"'"(2,]R) is a pair (G, C) G 

5L(2,M) X (0,oo), such that B = C-n(G), where H : 5L(2,M) SL{2,R) is the covering 

projection. Therefore all properties of lifts in SL{2, M) and their translation numbers are 
shared with the lifts Gb, B G GL+(2, M). 

Now we are able to show that if {xi,yi)i^z is a {p,q)-tjpe orbit of a positive twist 
map, F, then some matrices D(^Xi,y,)F G SL{2,M.), and hence their GL^(2, M)-conjugate 
matrices that represent the derivative in (x, x')-coordinate, cannot have a lift of trans- 
lation number in (0, 1/2), as it was supposed in [23]. Actually, if for some j, -D(a;^ 

has the trace in (—2,2), its basic lift, G, belongs to the subset S^i C SL{2,'R), i.e. 
t{G) G (—1/2,0), and any other lift G + j, j G Z, has the translation number in the 
interval (—1/2 + j, j). Thus no lift G + j can have r G (0, 1/2). 

Modifying correspondingly the Mather proof in [23], that is taking as lift for each 
factor in ( |T71) the basic lift, one gets the following result: 

Theorem 4.1 //(x„,|/„),n G Z zs a {p,q)-type orbit of the positive twist map F, Hg is 
the Hessian of the action Wpg atii.= (x„), / the Morse index of {xn), and r is the twist 
number of the corresponding f -orbit, then I and r are related as follows: 

i) If I is even, then t = —1/2; 

ii) If I is odd, then \—I/2~\ — 1 < r < [— //2], with equality if and only if ^ is a 
degenerate critical point. \a\ denotes the least integer number greater or equal to a G M. 

From Theorem 14. H and the general properties of lifts in SL{2, M) we have the following: 




B G GL+(2, M) ^ B, ^/det{B) G SL{2, R) x (0, oo) 

\^det{B) I 
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Corollary 4.1 i) The twist number of a minimizing orbit is zero. 

ii) Along an elliptic mini-maximizing q-periodic orbit, the maps f^, k = l,q have 
positive twist property. 

Proof: ii) Actually, if A; = 2, g violated the twist property, then the k- basic lift, 
Gk, associated to the orbit would have a displacement, such that ^^^(1/4) < —1/2, thus 
contradicting, in view of Lemma that the twist number is r G (—1/2, 0). 

Naturally, one raises the following questions: under what conditions a periodic orbit of 
twist number less than —1/2 can exist? Do such orbits so rarely appear that no reference 
to their existence in the dynamics of most studied twist maps can be found in literature? 

The Proposition below points out that such orbits typically appear at a period doubling 
bifurcation (details on this type of bifurcation can be found in [21], |25j): 

Proposition 4.1 Let f^, e > 0, be a positive twist map, which is a perturbation of an inte- 
grable map, /q. If at some threshold e' , a period doubling bifurcation of a mini-maximizing 
{p,q) -periodic orbit occurs, with transition from elliptic to inverse hyperbolic orbit, then 
the new born 2q-periodic orbit is elliptic, having the twist number within the interval 
(— 1, —1/2). // an inverse period doubling bifurcation occurs, i.e. with transition from an 
inverse hyperbolic to an elliptic q -periodic orbit, then at the threshold of bifurcation a new 
regular hyperbolic 2q-periodic orbit emerges, and its twist number is -1. 

Proof: We prove only the first part, the proof of the second one being similar. Let 
e I— J- Ze, e G (ei, £2), be the continuous path in the annulus described by the (p, g)-periodic 
orbit that at e' G (£1,62) undergoes a period doubling bifurcation. One associates to this 

path the continuous path, e ^ Gq{e), in the covering space, S'L(2,]R), described by the 
corresponding g-basic lift. The latter one starts in the region S-i (Fig. H]), intersects 
transversally the parabolic cone 'P2fc^_^i, and enters 'H'_i. 




SL(2,E) 



Figure 4: Illustration of the paths in the universal covering space SL{2, M), related to the 
period doubling bifurcation of a g-periodic orbit. 
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On the other hand, each is also a fixed point of /^^'', and thus the associated 2g- 
basic hft is G^(e). The twist number along the path e i— )■ G^(e) is double the twist number 
along the path e H- Gg(e), because r(G'^(e)) = 2r(G'g(e))), i.e. r(G^(e)) G (-1,0), for 
e G (ei,e')5 it gets —1 at e', and remains —1 within the interval (e',e2). Thus Gq(e) 
describes a path from E^2 to 7^-2, that intersects transversally the parabolic cone, ^2^^- 
From the point of intersection to V}2^ ^ a new path emerges, namely the path described 
by the 2g-basic lift K2g{e), associated to the new born 2g-periodic orbit, p^, e ^ {^',^2)- 

Because of the continuity of the maps e G [e',e2) — )■ — )■ K2q{e), and r : S'L(2,]R) — )• M, 
the new path enters the region 82- Thus the new born elliptic 2g-periodic orbit has the 
twist number within the interval (—1, —1/2). 

Because the 2g-basic lift associated to the new born elliptic 2g-periodic orbit has the 
range of its displacement map, [m{K2q), M{K2q)] C (—1, —1/2), each tangent vector at a 
point of that orbit is rotated clockwise in one period, with an angle between tt and 2tt. 

By Theorem 14.11 the twist number of a periodic orbit is well defined by an even Morse 
index of the corresponding critical sequence, but there is an ambiguity in deciding the 
demi-unit interval (an interval of length 1/2) that contains the twist number, when the 
Morse index is odd. For instance if the Morse index is 1, the twist number can belong 
either to [-1/2,0) or to (-1,-1/2). 

Because the translation number is not a group homomorphism (see E]), we cannot 
deduce the twist number of a {p, g)-periodic orbit, from the translation numbers of the 

basic lifts, D^f, . . . , 1^/9-1(2)/, associated to that orbit. 

If the Morse index of a {p, g)-critical sequence is an odd number, and the corresponding 
orbit is elliptic, we can use a naive method to get the demi-unit interval (i.e. an interval of 
the form {—{k + 1)/2, —k/2), k eN) containing its twist number. Such a method is based 
on the property that if z, z' are any two points on the same (elliptic in our case) {p, q)- 
periodic orbit, then D^f^ and D^'f^, are 5*17(2, ]R)-conjugated. Hence the corresponding 

g-basic lifts, Gg, Hg, are 5'L(2, M)-conjugated and by Proposition IB. H their displacement 
maps have the range included in the same demi-unit interval (— (fc + l)/2, —k/2), for some 
integer k >0. 

In order to get the actual value of k, we can choose an arbitrary point z on the elliptic 
g-periodic orbit, and compute a (g + l)-length segment of tangent trajectory, {zi,Wi), 
Wj+i = Dz-f{wi), i = 0,g, with zq = z and wq = (1,0)"^ G T^A. Let q be the x- 
coordinates of the vectors Wi, i = 0,q. If /c changes of sign in the sequence (q), i = 0,q, 
are recorded, then the twist number belongs to the demi-unit interval, (— (/c + l)/2, —k/2). 

This method works well for low values of the period q. In this way the reader can check 
easily the result in Proposition 14.11 in the case of standard map (Eq j28|) corresponding to 
e = 2, respectively e = 4, for which period doubling bifurcation of the 2-periodic orbit 
(0,0.5), respectively of the fixed point (0,0), occurs. 

The naive method is not suitable for hyperbolic periodic orbits, because in this case 
m{Gg) and M[Gq) belong to distinct demi-unit intervals. 

In the sequel, we present a general method to deduce the smallest interval X of ends 
in Z/2 that contains the interval [m{Gg), M{Gg)] of the twist number of a (p, g)-periodic 
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orbit. The method is based on the spectral properties of the Hessian matrix, Hq, and of 
an associated symmetric matrix. 

The Hessian matrix Hg of the Wpq-action associated to a {p, g)-periodic orbit, q > 3, 
is a Jacobi periodic matrix (i.e. a symmetric tridiagonal matrix with non-null entries in 
the upper right, and left lower corners, and the next-to-diagonal entries have the same 
sign). For q = 2, H2 is simply a symmetric matrix. 








ai 









/3i 








Oiq-2 Pq-2 
Pq-2 Olq-l 



g > 3, H2 



ao /3o + Pi 
/3o + /3i ai 



(19) 



where G M, /3j < 0, z = 0, g — 1, are defined in f[T8|) . Such a Jacobi matrix is called 
periodic because the sequences (a^), (/?«), i E Z, are g-periodic. 

The spectral properties of periodic Jacobi matrices were studied extensively in the late 
seventies. In his approach to the twist number, Angenent p] exploited a result by van 
Moerbeke [27] . The partial characterization of the spectrum of periodic Jacobi matrices 
given by van Moerbeke is complemented in [12], [7], and a thoroughly presentation of its 
properties, in connection to discrete Hill equation, can be found in [28|. 

Exploiting the properties of periodic Jacobi matrices we are able to visualize the results 
in Theorem 14.11 and moreover to enhance it. We give a method to deduce the smallest 
interval X of ends in Zi/2, containing the twist number of a {p, g)-periodic orbit, from the 
position of the real number zero, with respect to an interval whose ends are eigenvalues 
of the Hessian Hg, and/or eigenvalues of the symmetric matrix, H~, which is a slight 
modification of H„: 



/3o 




-Pq-l 





Pi ... 

Pq-3 Olq-2 

. . . Pq-2 



-Pq-l 




Pq-2 
Otq-1 



,q>3, 



«o Po - Pi 
Po - Pi ai 



(20) 



The Theorem 14. II can be enhanced to: 

Theorem 4.2 // {xn,yn),n E Z is a {p,q)-type orbit of the positive twist map F, Hq 
is the Hessian of the action Wpq at :s. = (x„), /, the Morse index of {xn), t, the twist 
number of the corresponding f -orbit, and I' is the number of negative eigenvalues of the 
symmetric matrix H~ , associated to Hq, then /, /' and r are related as follows: 

i) If I is even, then r = —1/2; 

a) If I is odd and Hq has at least one zero eigenvalue, then r = — (/ + l)/2; 
Hi) If I is odd and Hq has no zero eigenvalue, then: 
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a) If I' = 1-1, then re (-//2, -(/ - l)/2); 

b) Ifr = I, then t = -1/2; 

c) If r = 1 + 1, then re (-(/ + l)/2, -7/2); 

Proof: 

In order to relate the spectrum of the Hessian Hq to the trace of the derivative DzF'^, 
at a point z = {x,y) of a {p,q)-tjpe orbit, one associates to the Jacobi periodic matrix, 
J = Hg, the discrete Hill equation [28], = A^, A G M, where ^ = (^j), is a g-periodic 
sequence in M. The coordinate-wise version of this equation is: 



+ (ai - + = 0, AG 



(21) 



Among the solutions ,^(A) = (^i(A)) of this second-order difference equation, of par- 
ticular interest is the Floquet solution, i.e. the solution , that satisfies the condition 
^i+q{X) = ?7(A)^.j(A), with ^o(A),G(A) given, and //(A) an eigenvalue of the monodromy 
matrix: 



M(A) 







1 

X — Oq 





'(32 



1 

A - ^2 



1 

A) A - Qi 



(22) 



The matrix M(A) belongs to SL{2,R), while its factors to GL+(2,M). We note that 
for A = 0, M(0) is the derivative D^F'^, expressed in (x, x')-coordinates (see (fT7|) ). 
The trace of the matrix M(A) is a polynomial of order q, with real coefficients. 



tr(M(A)) = (/3o/3i---/3,-i, 



A« + 



(23) 



called the Hill discriminant |28] . 

The characteristic polynomial of a periodic Jacobi matrix, J, defining the discrete Hill 
equation is related to the Hill discriminant by the so called Hill formula [I2], [28] : 



det( J - XI) = {-iyf3of3, . . . /3,_i(tr(M(A)) - 2) 



(24) 



The history of Hill formula and generalization to multidimensional discrete and con- 
tinuous Lagrangian systems can be found in |6j. 

From (|24l) it follows that the eigenvalues, Aq, Ai, . . . , Ag_i, of the Hessian, Hq, are the 
roots of the equation tr(M(A)) = 2. The eigenvalues can be ordered in the following way 
[28]: 

A2i<A2i+i<A2^+2,V2>0 (25) 
On the other hand the roots of the equation tr(M(A)) = —2 can be ordered [28] as: 



and the two series are interlaced: 

Ao < Aq < X[ < Ai < As < A2 < 



(26) 



(27) 
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Being roots of the equations tr(M(A)) = 2, tr(M(A)) = -2, with M(A) G SL{2,R), Aai+i, 
respectively Agj, i > 1, can have at most double multiplicity. Thus each root occurs in 
(!27|) as many times as its multiplicity. 

For practical purposes it is important to note that the roots \[, i = 0, q — 1, are the 
eigenvalues of the symmetric matrix H~, defined in (1201) . associated to the Hessian Hq. 
Indeed, by |7] the characteristic polynomial of H~ is for any q > 2: 

det(if- - XQ = det{Hq - XQ + ■ ■ ■ 

Comparing to the expression of the Hill discriminant fl23l) one gets, that tr(M(A)) + 2 = 
iff det{H--XIg) = 0. 

For the reader convenience we illustrate in FigjS] the graph of a Hill discriminant 
corresponding to an even integer, q. 

Because the leading coefficient of the polynomial tr(M(A)) is (/3o/9i ■ ■ ■ Pg-i)~^, it fol- 
lows that in both cases, q even and odd, limA-j.-oo tr(M(A)) = +oo. Hence the dynamical 
type of the monodromy matrix (Figl5]) is as follows: 

a) If A < Ao or A G (A2i+i, A2j+2), i > with A2j+i ^ A2i+2, then M(A) is regular 
hyperbolic. 

b) For A G (A2i, AsJ or A G (A'aj+i, A2i+i), i > 0, M(A) is elliptic. 

c) For A G (A'ai, X'^i+i), X'^i ^ A'a^+i, i > 0, M(A) is inverse hyperbolic. 

Now the location of the real number zero with respect to the interlaced series of 
eigenvalues (1271) gives information on the Morse index of Hq, the signature of the matrix 
H~, and the dynamical type (regular hyperbolic, elliptic or inverse hyperbolic) of the 
{p, g)-periodic orbit we are studying. 

If the Morse index / of the Hessian Hq is odd, and Hq has no zero eigenvalue, i.e. 
G (A2i,A2i+i), for some i > 0, then by Theorem 14.11 the twist number belongs to an 
interval of length 1, \-I/2] - 1 < r < \-I/2]. It follows that the g-basic lift Gg 
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associated to our orbit can belong to £-2i+i,H'_2i+ii or to £-21-1 C 5*^(2, M), and the 
twist number can belong to different subintervals of the above 1-length interval. In order 
to deduce the right demi-unit interval to which the twist number belongs we analyze 
the position of with respect to the eigenvalues Agj, Agj+i, and deduce the index I' of the 

symmetric matrix . Due to the continuity of the translation number r : SL{2, M) — > M, 
the case iii) from the Theorem follows directly. 
Similarly one can prove the case ii). 

Corollary 4.2 The smallest twist number a {p, q) -periodic orbit of a positive twist APM 
can have is —q/2. 

This property was also stated in [2] but no arguments were given. 

In numerical investigations of the twist number, if the Lagrangian generating function 
of the twist map is defined only locally (as in the case of tokamap |4|), then one can 
reduce the computations and improve accuracy, evaluating the Hessian Hg associated to a 
periodic orbit, not in (x, x')-coordinates, but in the original symplectic (x, y)- coordinates. 
In this case the entries of Hq are as in ( |TT|) . 

Thus, a simple algorithm that in the first step calculates the eigenvalues of the Hessian 
( ITTj) . order them increasingly, and then tests the cases described in Theorem 14.21 (com- 
puting also the eigenvalues of H~ , if necessary, and ordering them increasingly), returns 
the value of the twist number or the demi-unit interval containing the twist number. 

We note that the 2g-periodic orbits born through a period-doubling bifurcation are 
ordered orbits (Birkhoff orbits). In the next section we give examples of badly ordered 
(non-Birkhoff ) {p, g)-periodic orbits that are born through a saddle-center bifurcation 
and have the absolute twist number greater than 1/2. We also show that the so called 
two-component strong folding region maps can exhibit ordered periodic orbits of large 
absolute twist number. 



5 Twist number and order properties of periodic or- 
bits 

Computer experiments reveal that plenty of pairs of periodic orbits of absolute twist 
number, |r| > 1/2, can be found in the folding region of the standard map corresponding 
to the perturbation parameter e > 1. 

The standard map we consider here is defined by the lift F : (x, y) (-)■ (x', y'): 

X + y' 

y - ^sm{27Tx), e G M, ^'^^^ 

Example 1. The standard map corresponding to the perturbation parameter e = 
1.4578 has the minimizing (8, 3)-periodic orbit represented by (0.5,0.693821664066481), 
and the mini-maximizing one by (0,0.586319735666097) (Fig. [6l left panel). 
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Figure 6: Four (8, 13)-periodic orbits of the standard map corresponding to e = 1.4578. 
At left is illustrated the pair of orbits given by Aubry-Mather theory, and at right a pair 
of orbits of twist number —1, respectively rG(— 1,-1/2). 



The point z = (0.5, 0.7237177352891645) is also a (8, 13)-periodic point. Applying the 
ubiquitous method used in the study of twist area preserving maps, namely, computing 
only the residue, R = 0.44794945246558, we can mistakenly conclude that it is a mini- 
maximizing periodic orbit. Computing further the eigenvalues of the Hessian Hi^ we get 
the Morse index 1 = 1, and analogously we are tempted to conclude that it is mini- 
maximizing (the index equal to one is only a necessary condition to be mini-maximizing, 
not a sufficient one). But the orbit is unordered and thus it cannot be mini-maximizing. 
The eigenvalues of the associated matrix H^^ reveal that its twist number belongs to 
(—1,-1/2). Hence this orbit is an elliptic periodic orbit with the property that each 
vector V G T^A is rotated clockwise about the orbit of z, more than vr radians in one 
period. 

The orbit of z was born through a saddle-center bifurcation at a slightly lower value of 
e, and its pair is the regular hyperbolic orbit of the point z' = (0.5,0.723372427461496). 
z' has the twist number r = — 1, and its orbit is also unordered. They are illustrated in 
Fig. ini right panel. 

One can observe that projecting each of the four orbits onto the first factor, 
identified with [—0.5,0.5), and denoting by X£, Xr the abscissas of the minimizing-orbit, 
nearest to {xe at left and Xr at right of 0), then within the interval {xe, Xr) each projected 
non-Birkhoff orbit has two points. 

Example 2. The point z = (-0.387429398213243, 0.225141203573513) of the stan- 
dard map corresponding to £ = 3 has an inverse hyperbolic (1, 7)-periodic orbit of trans- 
lation number r = -1.5, and its pair z' = (0.389952663038891, 0.536877082952284) has a 
regular hyperbolic (1, 7)-periodic orbit of twist number r = — 1. 

For the same parameter the Aubry-Mather orbits of type (1,7) are the following: 
(0,0.390375216334041475), mini-maximizing, and (-0.5,0.008258333079127612), mini- 
mizing. 

Considering again the interval {xg, Xr) as in the previous example, we note that within 
this interval each orbit of large absolute twist number has three points. 
In FigJ7]are illustrated all four orbits of type (1, 7). 

In an attempt to asses numerically the efficiency of the Boyland converse KAM crite- 
rion [8], Leage and MacKay [19] looked for badly ordered periodic orbits of the standard 
map. They detected such orbits for large perturbations parameters. 
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Figure 7: Four (1, 7)-periodic orbits of the standard map corresponding to e = 3. At 
left is illustrated the pair of orbits given by Aubry-Mather theory, and at right a pair of 
orbits having the twist number —1, respectively —1.5. 

We note that the badly ordered orbits reported in [19] have also non-zero twist number 
and |r| > 1/2. They are orbits born through a saddle-center bifurcation, at a parameter 
value greater than the threshold of the period doubling bifurcation of the mini-maximizing 
periodic orbit of the same type (p, q). 

Example 3. The point z = (0.181826060531142,0.681826060531142) found in [I9] 
for e = 7.221365 is an elliptic point of type (1,3), twist number r G (—1, —1/2), and its 
pair z' = (0.181803685309653,0.681803685309654) has a regular hyperbohc (l,3)-orbit, 
and r = —1. 

The above examples can suggest that (p, g)-periodic orbits of large absolute twist 
number are unordered. A natural question is whether a positive twist map can also 
exhibit ordered {p, g)-periodic orbit [p, q relative prime integers) of twist number less 
than —1/2. The following example shows that the answer is positive. 

Example 4. We consider a multi-harmonic standard map defined by the following 

lift: 

x' = x + y + P{x) , . 

y' = y + Pix), ^^^^ 

where P is the trigonometric polynomial, P{x) = rji sin(27ra;) + 772 sin(47ra;) + 773 sin(67rx), 
corresponding to the point rj = (0.18, —0.42, —0.11) in the parameter space. This map is 
reversible and has four (1, 2)-periodic orbits (FigJHl third line, right) : 

a) (0.220481551875563,0.440963103751125) is a symmetric minimizing orbit. 

b) (0, 0.5) is a regular hyperbolic point having the twist number equal to —1 (the associ- 
ated Hessian has negative eigenvalues Ai = -5.48881762834532, A2 = -1.06693368771638, 
and thus the corresponding (1, 2)-sequence is a sequence of maximum action). 

c) zi = (0.099868507451075,0.699123941758688), 
22 = (0.900131492548925,0.300876058241312) 

are two inverse hyperbolic points having distinct orbits, of twist number —1/2. 
All four (1, 2)- periodic orbits are cyclically ordered. 

We are wondering why in this case an orbit of large absolute twist is ordered, while 
in the case of standard map such orbits are badly ordered. In order to give an answer we 



22 



connect this map to an integrable one, i.e. we consider the maps F^, e G [0, 1] defined by: 

x' = x + y + eP{x) , . 

y' = y + ePix), ^^^^ 

and we analyze the dynamical type and twist number of (1, 2)-periodic orbits, as e in- 
creases from to 1. 

The map F^, close to the integrable twist map, has two (1, 2)-periodic orbits: a min- 
imizing one that persists down to e = 1, and the orbit of the point (0,0.5) which starts 
as a mini-maximizing orbit of elliptic type (FigJHl first line, left). The latter one passes 
through a sequence of bifurcations (FigJHI) that lead to the following change of its twist 
number (see FigjO]): 

elliptic mini-maximizing — )■ inverse hyperbolic ofr = — 1/2— > elliptic of twist number 
in (—1,-1/2) — 7- regular hyperbolic of r = — 1. 

The first two bifurcations are period doubling, while the last one is a Rimmer-type 
bifurcation that leads to the creation of the two distinct asymmetric periodic orbits, 0{zi), 

0{Z2). 

Hence the maximizing orbit (0, 0.5) of the map Fi was born as an ordered mini- 
maximizing orbit, and through a sequence of bifurcations it changed the twist number 
as it is illustrated in FigJH but not the order property, because the period doubling, and 
Rimmer bifurcation of ordered orbits leads to ordered orbits. 

We also note that while the (1, 2) -minimizing orbit does not undergo any bifurca- 
tion as the perturbation parameter increases, the (l,2)-mini-maximizing orbit bifurcates, 
changing also the twist number value after each period-doubling bifurcation. Moreover at 
the Rimmer bifurcation threshold, the symmetric mini-maximizing orbit is continued by 
a maximizing one, and two asymmetric mini-maximizing orbits are born. Thus at e = 1 
there exist 3 orbits given by Aubry-Mather theory: a symmetric minimizing orbit, and 
two asymmetric mini-maximizing orbits. 

Remark 5.1 From the analysis of the map\3^it follows that the mini-maximizing {p, q)- 
periodic orbit of a twist map can have the twist number r either in the interval [—1/2, 0) 
or in the interval (—1,-1/2), while the twist number of the corresponding minimizing 
orbit is constant and equal to 0. 

In the sequel we show that the unexpected behaviour of the map (150]) is due to the 
fact that it is a TSF map. 

First we show that if a twist map is either USF or TSF we can predict the bifurcations 
that the orbits intersecting Fq can undergo, as the perturbation parameter varies. 

Proposition 5.1 Let F^ be an USF map for e > e* , and F* the symmetry line that 
get included in the strong folding region for e > e* . If after a slight perturbation of the 
integrable map Fq, a {p,q) -periodic orbit, q > I, has a point {xo{e) , yQ^e)) on T* , then the 
orbit can be either mini-maximizing or minimizing. Increasing e from 0, the orbit can 
undergo the following bifurcations: 
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Figure 8: The phase portrait of the map F^, Eq. (130|) . illustrating the sequence of bifur- 
cations that the orbit of the point (0, 0.5) undergoes as e increases from to 1. First line: 
at left, e = 0.028, (0, 0.5) has an elliptic mini-maximizing orbit; at right, e = 0.4, slightly 
above a period doubling bifurcation. Second hne: at left, e = 0.46, slightly above the sec- 
ond period doubling bifurcation (from r = —1/2 to r G (—1, —1/2)); at right, e = 0.65, 
(0,0.5) has an elliptic orbit of twist number r G (—1, —1/2). Third line: at left e = 0.8, 
after a Rimmer bifurcation (from r G (—1, —1/2) to r = —1); at right, e = 1, (0, 0.5) has 
regular hyperbolic orbit of twist number r = — 1, and the periodic orbits born through 
Rimmer bifurcation (the orbits marked with -|-, respectively x, are inverse hyperbolic of 
r = -l/2). 




0.25 0.5 0.75 
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Figure 9: The twist number of the (1, 2)-periodic orbits map (Eq. [301) . as a function of 
e. The dashed line is the twist number of the minimizing orbit, the continuous line is the 
twist number of the point (0,0.5), and the dotted one is the twist number for the two 
asymmetric periodic orbits, born through a Rimmer bifurcation. 
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a) If {xQ{e),yQ{e)) is mini-maximizing of elliptic type it undergoes necessarily a period 
doubling bifurcation at some < ed < e* . 

b) If {xo{e),yo{e)) is minimizing, then at a threshold e.^ it gets mini-maximizing through 
a Rimmer-type bifurcation, and at some > e,- it undergoes a period doubling bifurcation, 
and gets inverse hyperbolic of twist number r = —1/2. 

Proof. 

At e = S{x,y) = 2, for every point {x,y). Let Hg be the Hessian matrix associ- 
ated to the (p, g)-orbit, having a point on T*, H~ the symmetric matrix (l20l) . and Aq, 
respectively Ag, the smallest eigenvalues of these matrices. From Remark 13.11 it follows 
that the function 6{xq, yo; e) decreases from 2, as e increases. By Lemma [XT] we have that 
Ao, Aq < S{xo,yo) < 2, and by Gershgorin Theorem Aq belongs necessarily to the interval 
[S{xQ,yo]e) — 2,6{xQ,yQ]e)]. Hence for e very close to zero, Aq can be either negative or 
positive, i.e. after a slight perturbation of the integrable map the considered periodic 
orbit can be mini-maximizing of elliptic type or minimizing. 

a) If the orbit is mini-maximizing, then since the function 6{xo, yo; e) is decreasing, 
the eigenvalue Ao(e) is pushed to the left on the real axis, and at a threshold < e* it 
crosses 0, i.e. the twist number crosses —1/2 and a period doubling bifurcation occurs. 

b) If the orbit is minimizing, using similar arguments it follows that necessarily the 
stated sequence of bifurcations must occur. 

Corollary 5.1 An USF map exhibits a total or a partial dominant symmetry line, that is 
either for each perturbation parameter or for e greater than a threshold, e* , any periodic 
orbit having a point on Fq has the twist number less than zero. 

It is conjectured that the standard map has Fq as a total dominant symmetry line 
(as far as we know, no proof was given for this numerical observation). The map (Eq JT^ 
corresponding to a = —3 has Fq as a partial dominant symmetry line. 

In the case of a TSF map a similar behaviour have the orbits intersecting Fq, but 
moreover a second period doubling bifurcation from r = — 1/2 to r G (—1, —1/2) is also 
possible. 

Let us analyze the (1, 2)-periodic orbit of the twist map fl30|) . intersecting Fq. 

After a slight perturbation of the integrable map the mini-maximizing (1, 2)-periodic 
orbit has its points (xo,2/o), (3^1,1/1) on the symmetry lines Fo,Fg. The restriction of the 
1-cone function 5 to each of these symmetry lines is a decreasing function with respect 
to e, namely it decreases from the value 2 corresponding to e = 0, to negative values 
corresponding to e for which Conei(e) has two connected components. It follows that as 
e increases, the smallest eigenvalue, Aq < min{(5(xo, |/o), |/i)}, of H2 is pushed to the 
left on the real well as the smallest eigenvalue Aq = min{(5(xo, 1/0), ^(a;!, |/i)} of 

the matrix H2 (see [2UI) . Hence a sequence of bifurcations of the ordered (1, 2)-periodic 
orbit occurs leading to a continuous decrease of the twist number. Beyond an e both 
points of the orbit are within the strong folding region. In this case the Hessian H2 and 
its associated matrix H2 have negative diagonal entries, and as a consequence H2 has 
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Figure 10: The profile of S function associated to the standard map (left), respectively a 
TSF map (right) and the value of 6 at mini-maximizing points of the (1, 2)-orbit (marked 
by o), and minimizing points (marked by *). 



both eigenvalues negative. Thus the (1, 2)-periodic orbits intersecting Fix(R) is in this 
case an ordered orbit, of twist number r < —1/2. 

Let us show that in the case of an USF map the second period doubling bifurcation 
of a mini-maximizing (1, 2)-periodic orbit is not possible (i.e. transition from r = —1/2 
to r G (-1,-1/2)). 

By Proposition 15.11 The points of the mini-maximizing (l,2)-orbit of an USF map 
are forced (eventually after a threshold) to belong to the symmetry lines FcFq. The 
points, (xq, 1/0)5 i^i^y'i)} of the minimizing periodic orbit belong to Fi, respectively F'^, 
and i?(xo),?/o) = {x'i,y[). Taking into account that for standard-like maps the 1-cone 
function is /^-invariant, i.e. 

SiR{x,y)) = S{x,y), 

for any point {x, y) on a symmetric periodic orbit, we have the following Hessian matrices 
(fT9|) if2, associated to the two orbits, and the corresponding symmetric matrices (120|) 
H2, H'2 : 



H, 



S{xo,yo) 
-2 



-2 

Sixi,yi) 



6{xo,yo) 
S{xi,yi] 



(31) 



H' 



^ix'Q,y'Q) 
-2 



-2 
^ix'o,y'Q) 



H'- 



^ix'Q,y'Q) 






^ix'o,y'Q) 



(32) 



In Fig{TD] is illustrated the profile of the 6 function for the standard map, as an USF 
map, and for the TSF map fl30p . as well as the values of S along the mini-maximizing, 
respectively minimizing (1, 2)-periodic points of these maps. 

A simple algebraic computation reveals that the different relative position of diagonal 
entries of H'2 with respect to those of H2 , in the two cases, leads to distinct relative 
position of the corresponding eigenvalues (Fig JTTl) . In the case of the USF map (FigJTTl 
left) the period doubling bifurcation from r = —1/2 to r G (-1,-1/2) is obstructed 
because the double eigenvalue of H'2 is positive, and less than the second eigenvalue, 
\,iH2),ofH2. 
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Figure 11: The graphs of the Hill discriminants associated to (1, 2)-periodic orbits of the 
standard map (as a USF map), corresponding to e = 1.3 (left) respectively of the TSF 
map f pOj) . corresponding to e = 0.8 (right). These maps have the profile of 6 function 
illustrated in Fig. [TOl 

In the case of the TSF map (FigJTH right) as the perturbation parameter increases, 
the double eigenvalues of H'2 increases, while Xo{H2), \i{H2) decrease and successively 
they get zero, without any obstruction. 

We conjecture that in the case of an USF map no mini-maximizing {p, q')-type orbit of 
twist number —1/2 can bifurcate to a an orbit of twist number rG(— 1,-1/2), because of 
a similar obstruction. In this case the only periodic orbits of large absolute twist number 
are born through a saddle-center bifurcation. 

After the discussion concerning the residue and the index of the Hessian matrix asso- 
ciated to periodic orbits given in Example 1 we are tempted to state that among residue, 
the Morse index, and the twist number of a (p, g)-periodic orbit, the twist number is the 
only numerical quantity that leads to a correct decision whether or not that orbit is one 
given by Aubry-Mather theory. However this is not the case, because as the following 
example illustrates, there can exist badly ordered periodic orbits of twist number zero. 

Example 5. The point z = (0.5, 0.323020669964897189) of the standard map corre- 
sponding to e = 2.87 is a (3, 5)-regular hyperbolic point, because its Hessian has the 
eigenvalues: 

0.04121176201302, 0.53521014626564, 3.88911846604733, 
5.36539881978110, 5.89051890157580 

Its twist number is zero, but it is not minimizing because it is badly ordered. 

The Hessian associated to its pair, z' = (0.4941763525132936290, 0.3344329777751245200), 
has the eigenvalues: 

-0.07386797172616, 0.56782907421433, 3.88372973742859, 
5.35939921395639, 5.88621274711297 

The orbit of z' is also badly ordered. The minimizing and mini-maximizing periodic orbits 
of the type (3, 5) are respectively: 

(0.5, 0.792447770847267734), (0, 0.567936349180897305) 
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The action along the orbits of these four periodic points of type (3, 5), in the order of 
their mention, is respectively: 

0.891078196747568, 0.891083047652180, 0.804181852816563, 0.831286899496951 

6 Conclusions 

Deciphering the structure of the universal covering space of the group S'L(2,M), and 
exploiting the properties of the translation number of circle maps, we developed a general 
framework for the study of the twist number of periodic orbits of twist maps. 

We developed an algorithm to compute the value of the twist number or an interval 
containing this value. Applying this algorithm we illustrated that the classical method 
used in dynamical systems theory, to classify periodic orbits of a twist map only by their 
residue or the Morse index of the associated Hessian matrix can lead to errors. 

We proved that ordered periodic orbits of large absolute twist number are born through 
a period doubling bifurcation. We also gave examples of periodic orbits born through a 
saddle-center bifurcation that have large absolute twist number and are unordered. 

In order to explain the variation of the twist number in a subclass of standard-like 
maps we introduced a new tool of study of the dynamical behaviour of these maps, namely 
the so called 1-cone function. The location of minima of this function, as well as its values 
along a periodic orbit appears to give valuable information on that orbit. 

The examples given in this paper illustrate that a correct characterization of a periodic 
orbit of a twist map needs the analysis of the residue, the twist number (i.e. the sign of 
eigenvalues of the Hessian matrix, H, and of its companion matrix H~), as well as the 
order properties of that orbit. 

The twist number allows a better characterization of the periodic orbits given by 
Aubry-Mather theory: a minimizing orbit is ordered and has the twist number equal to 
zero, while a mini-maximizing orbit is also ordered, but its twist number can belong either 
to the interval [—1/2,0) or to the interval (—1, —1/2). 

The study initiated in this paper allows to investigate the rotational component of 
chaos, that along with the Lyapunov component can give a deeper insight into the dy- 
namics of twist area preserving maps. 

The present study is not exhaustive and implications of the existence of orbits of large 
twist numbers are to be studied further. We mention here only few problems that deserve 
to be investigated: 

a) What particular properties of a twist area preserving map lead to the birth of badly 
ordered periodic orbits of large absolute twist number. 

b) How invariant manifolds of hyperbolic periodic orbits of absolute twist number 
greater or equal to 1 intersect with the invariant manifolds of minimizing periodic orbits 
and how such intersections influence the transport in the phase space or the dynamics in 
a region of instability. 
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A Translation number of circle homeomorphisms in- 
duced by matrices in SL(2^ M) 

A.l Translation number of an orientation preserving circle home- 
omorphism 

Our approach of the twist number of a g-periodic orbit of a twist APM exploits some 
properties of orientation preserving homeomorphisms (OPHs) of the unit circle, and of 
the translation number of their orbits. For details concerning dynamics of OPHs of the 
circle we refer to [HI [T71 [13] . 

In the sequel we consider §^ either as the unit circle, C, in or as the quotient group 
M/Z. The two definitions are equivalent through the identification: 

e + Ze R/Z ^ ue = (cos(27r^), sin(27r^)) G C C (33) 

We denote by tt : M — S"*^, the covering projection, tt{x) = x (mod 1). A lift of an 
OPH : — )■ is a homeomorphism, G : R — M, such that tt o G = g o tt. The 
orientation property of g ensures that G{x + 1) = G{x) + 1, V x G M. Thus the function 
= G — id]8 is 1-periodic and can be written as = ip o tt , where : §^ — t- R. is called 
the displacement map of g associated to G [13]. The value of ■?/' at a point of measures 
the amount that point is displaced by g around the circle. 

In the sequel we especially use the function "if = G — id^, and call it the displacement 
map of G. 

We denote by m(G) = min^gK ^(a^), M{G) = maXj;gR^(x). 

Two lifts G,H of g differ by an integer, H = G + k, as well as their displacement 
maps, = "^G + k, k e Z. 

The group Homeo'''(§^) of OPHs of the circle has the universal covering group [2U] : 

Homeo^(§^) = {G : R — t- R | G strictly increasing continuous map, G{x+1) = G{x) + 1}, 

i.e. the elements of this groups are lifts of OPHs of S^. We denote by H the covering 
projection: 

n : Homeo+(§i) ^ Homeo+(§^), H(G) = g, 

where g is the unique homeomorphism in Homeo^(S^) such that gn = tiG. H is a group 
homomorphism. 

An important role in the qualitative description of the dynamics of OPHs of the circle 
has the Poincare translation number. 
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The translation number of a homeomorphism G G Homeo^(§^) is defined by: 

ri^G) = lim ^"^""^"^ 

n— i>oo n 

This number exists and is independent on x. 

For a self-contained presentation we list some properties of the homeomorphisms in 

Homeo'''(§^) and their translation numbers: 

Proposition A.l Let g : E>^ ^ be an OPH, and G a lift of g. Then the following 
properties hold: 

1. The translation number is invariant under the conjugacy in Homeo^ {S^) , that is 
t{H-^ oGo H-^) = t{G), for every G,H e iJomeo+(§i). 

2. The translation number r : Homeo^ (E>^) — t- M zs not a homomorphism of groups, 
that is, in general 

r(G2oG'i)^r(G2) + r(Gi), 

and its deviation from being a homomorphism is bounded, existing a positive number D 
such that: 

\t{G2 o Gi) - t{G2) - r(G'i) < Z^, V Gi, G2 G i/or?^(§i). (34) 

The least constant D with this property is 1 [29] . 

3. T^G") = nriG), for any positive integer n. 

4- The translation number of the inverse of a homeomorphism G G Homeo'^ {S>^) is 
r{G-') = -r{G). 

5. The difference between the max and min displacement of a lift G is less than 1, 
M{G)-m{G) < 1 [IE]. 

6. The translation number of a homeomorphism, G, belongs to the range of the dis- 
placement map, T G [m{G) , M (G)] . 



A. 2 The translation number of the circle homeomorphisms de- 
fined by matrices in SL{2, R) 

SL{2,M.) is the group of 2 x 2 real matrices of determinant equal to 1. It acts on the 
vector space M^. This action induces a natural action of the group S'L(2,R) on the 
unit circle, S^, i.e. there is an injective homomorphism : 5'L(2,R) — t- Homeo^(§^) 
between SL{2,M.) and the group of orientation preserving homeomorphisms of the circle 
C = {ue = (cos(27r^), sin(27r0) G R^^ ^ G M}: 

A G SL{2, R)^ gAE Homeo+(§^), (35) 

where gAiug) = - — and || ■ || is the Euclidean norm in R^. 

\\Aue\\ 
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Therefore the group S'L(2,M) is identified with a subgroup of Homeo'''(S^), and its 

universal covering group S'L(2,M) with a subgroup of Homeo^(Si). We denote by 11 : 

SL{2,M.) —7- SL{2,'R) the covering projection. 

Our aim is to characterize the translation number of the lifts G G SL(2, M). Its particu- 
lar properties in comparison to the translation number of a general circle homeomorphism 
derive from the following : 

Proposition A. 2 Any lift Ga of an OPH, gA,A E S'L(2,M), has the property 

Ga{x + 1/2) = GAix) + 1/2, y X eR (36) 

Proof: We denote by T^, a; G M, the shift homeomorphism of M, T^{x) = x + u. The 
center of the group SL{2,'R) is {±/2}. Because each matrix A commutes with —I2, the 
circle map, qa, commutes with the circle homeomorphism, h, defined by: 

h:9 + ZER/Z^e + l/2 + ZER/Z (37) 

The lifts of this map are the shifts T{2k+i)/2i k E Tj. Thus a lift Ga commutes with 
Ti/2, and this amounts to G{x + 1/2) = G{x) + 1/2, V x G M. 

Let Homeo^(S^) be the subgroup of orientation preserving homeomorphisms of the 
circle, that commute with the map h defined in (EZ]). Because SL{2, M) C Homeo^(§^), we 
deduce the particularities of the translation number of the homeomorphisms in Homeo^(§^) 

and those properties will be inherited by the translation numbers of lifts in SL{2, M). The 
key point in deriving these particularities is that Homeo^(S^) is a rescaled version of 

Homeo'^(S^). Indeed, let zu : Homeo'''(S^) — t- Homeo^(S^) be a map defined by: 

w{H){x) = ^H{2x),x E R (38) 

tu is a monomorphism of groups, and its image is Homeo^(§^). Hence each G E Homeo^(S^) 

is represented by a general homeomorphism H E Homeo'''(§^), H = zxj^^{G), H{x) = 
2G{x/2). 

Proposition A. 3 LetG be a homeomorphism in Homeo'j^ (Ei^) . Then its translation num- 
ber is half the translation number of the corresponding rescaled version H = w~^{G): 
r{G) = t{H)/2. 

Proof. From G{x) = -H{2x), we get G'^{x) = -iJ"(2a;), and the translation number 

is: 

riG) = i lim = 

^ ' 2n^oo n 2 

From the definition of tu, the Propositions IA.2 1 IA.3t and the properties of general circle 
homeomorphisms, we have the following properties of the displacement maps of lifts of 
circle maps defined by matrices in S'L(2,M), relevant to our approach to twist numbers: 
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Figure 12: The graph of displacement maps associated to two homeomorphisms Ga-,Ha-, 
Ha = w-\Ga) ■ 

Corollary A.l a) The displacement map of a lift Ga ^ S'L(2,]R) is a periodic function 
of period 1/2 (Fig\EM, left). 

h) If = Ga — id-R is the displacement map of a lift Ga G SL{2,'R), and $ is the 
displacement map of the lift Ha = ta''^{GA), then = -$(2x), V x G M (Fig. \A.^) . 

c ) The difference between the max and min displacement of a homeomorphism Ga G 
SLiX^) is: 

< M{Ga) - m{GA) < 1/2 (39) 

d) The restriction of the translation number r to the universal covering space of the 
group SL{2, M) has the property: 

\t{Ga o Gb) - t{Ga) - t{Gb) I < 1/2, V Ga, Gb e SLi^R) (40) 



B The structure of the universal covering space of 

5L(2,IR) 

We define a system of coordinates on the universal cover group S'L(2,M), that allows us 
to understand and visualize the topology of this group. 

We recall that the matrices in SL{2, M) are classified according to their trace. Namely, 
a matrix A G SL{2,M.) \ {±/2} is called regular hyperbolic, inverse hyperbolic, elliptic, 
positive parabohc, respectively negative parabolic if ti^A) > 2, ti^A) < —2, |tr(A)| < 2, 
tr(A) = 2, respectively tr^A) = —2. 
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Each matrix A G SL{2,M.) is 5'L(2, M)-conjugated with one of the normal forms: 



It is well known that two conjugated matrices have the same trace. The conversely is 
not always true. The matrices within the pairs, {R2-Ke, R-2n9), {P^'^, P^'~^), {.P~^'^, P~^'~^) 
for example, have the same trace, but they are are not S'L(2, ]R)-conjugated. 



is S'L(2, M)-conjugated to one of the above normal forms. 

a) An elliptic matrix, A e SL{2,M.), is S'L(2, M)-conjugated with a unique matrix 
i?27r6i- More precisely according to whether 6 — c>0or6 — c<0 (conditions that 
are equivalent respectively to 6 > or 6 < 0), it is conjugated to R2TTe, Q G (~l/2,0) 
or to R2-Kei d G (0,1/2) (this property can be checked by a straightforward matricial 
computation). A matrix with 6 = cannot be elliptic. 

Thus SL{2, M) contains two subsets of elliptic type matrices, one consisting in matrices 
that are SL{2, M)-conjugated to a clockwise rotation (or equivalently, they have the entry 
6 > 0), and another is formed by matrices conjugated to an anticlockwise rotation (their 
6-entry is negative). 

b) If A G SL{2, M) is regular hyperbolic, respectively inverse hyperbolic, then A is 
conjugated to an H^, respectively an iJ*-type matrix. 

c) If A G SL{2, M) is a positive (negative) parabolic matrix, then it is conjugated 
with the first or the second matrix of the pairs (-P^'^, P^'^^), (P^^'^, P^^'^^), according to 
whether 6 — c>0or6 — c<0. 

We can conclude that: 

Two hyperbolic matrices are SL{2,'R) -conjugated iff they have the same trace. Two 
elliptic, matrices, are SL{2, M.) -conjugated if and only if they have the same trace and the 
same sign of the entries b, and two parabolic matrices are SL{2,M.) -conjugated iff they 
have the same trace and the same sign of the differences b — c, associated to each matrix. 

Topologically, SL{2, M) is the interior of a three dimensional solid torus (IH], pag. 14) 
i.e. it is homeomorphic to D x S^, where D is an open disc centered to the origin in M^. 

Hence its universal covering group, SL{2,'R), is homeomorphic(as well as diffeomorphic) 

to DxM.. Obviously, SL{2, M) is also diffeomorphic to M^, and a few such diffeomorphisms 
are constructed in [9] . We define a diffeomorphism that fits best to our purpose of studying 





(41) 



In the sequel we point out under what conditions a matrix A 



a b 
a d 



G 5L(2,M) 
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the twist number of orbits in area preserving maps. More precisely we choose a system of 

coordinates, such that the visual representation of SL{2,'R) in this system looks like the 
Penrose diagram of S'L(2,M), regarded as an anti-de-Sitter space p. 

Any matrix A G 5*^(2, M) decomposes uniquely as A = RS (polar decomposition), 
where R = R2-Ke is a rotation matrix, and 5" is a symmetric positive definite matrix. Let 
D be the open disc of radius 1/2, centered at the origin in R^, and (p, a) polar coordinates 

within this disc. We define x : ^ x M ^ 5L(2, M), such that x(0, 0, 0) = id^ e SL{2, M), 
and the polar decomposition of the matrix A = n(x(p, a, 6)) be: 



A — RottqR 



sin (ttp) ^ 



1 + sin(7rp) 

' 1 + sin (vrp) 
1 — sin(7rp) 







R-2.a (42) 



One can prove that x is a diffeomorphism. Since the trace of the polar decomposition, 
RS, is trace(i?)trace(5')/2, we get: 

NN 2cos(27r0) , 

trace n x^, P, a = 43 

cos(7rp) 

In the sequel the trace of a homemorphism G G SL{2, M) is meant as the trace of the 
corresponding projection, Il{G), onto SL(2, M). According to their traces we also call the 

lifts, elliptic, hyperbolic or parabolic lifts. Denoting by tr : SL{2, M) — t- M the function 
trace, the connected components of the level sets Lt = tr^^{t}, t G M, are 2-dimensional 

surfaces and each such a component represents a conjugacy class in SL{2,M) [9]. 

From (H3|) it follows that tr does not depend on a, and thus each connected component 
of a level set, Lj, is a surface of revolution about the 9 axis. This property implies that 
the relative position of different conjugacy classes can be visualized in any planar section 
of the infinite solid cylinder, through a plane containing its symmetry axis (Fig JTS]) 

In order to analyze a section in the solid cylinder, we locate the particular lifts, that 

project to normal forms in SL{2,'R). We denote by T^^ G SL(2,M.), the shift homeomor- 
phism, defined by T^(x) = x + u. 11 (T^^) = -R27raj- 

The center of 5'L(2,M), i.e. the set, Z{SL{2,M.)), of the lifts that commutes with any 
other lift, consists in the set of shifts Tk/2, k E Z. n(Tfc/2) = ^h, according to as k is even 
or odd. n(x(p, a, /c)), /c G Z is a regular hyperbohc matrix, and n(x(p, a, {2k + l)/2)), 
A; G Z, is an inverse hyperbolic matrix. 

The components of the level set L2, and L„2 i-e. the set of positive, respectively 
negative parabolic lifts are cones in D x M with the vertex removed. 

Analyzing the regions inside the infinite open solid cylinder, defined by the two families 

. . , . , r 2cos(27r^) 

of cones, and the equations of the conjugacy classes surfaces, — — = t, correspond- 

cos(7rp) 

ing to t G (—2, 2), t > 2 and t < —2 we can conclude that: 
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SL(2,m) R J,) 

Figure 13: a) Visualization of the topological structure of the group 5*17(2, M) and the 

translation number function r : SL{2, M) — M; b) The foliation of the regions S-i, Hq, 'H'_i 
by conjugacy classes. 



1. The parabolic cones have the vertices at a shift Tfc/2, k G Z. More precisely, if k 
is even (odd), then is the vertex of a positive (negative) parabolic cone. Each set 
{Tk/2\ is a conjugacy class. 

2. Inside a region bounded by two patches of nearby cones with vertices at 7^/2, 
T(fc+i)/2, A: G Z, lie elliptic lifts. These regions are denoted by £k- For k even integer, the 
matrices Il.{£k) have 6 < 0, while for odd k they have 6 > 0. 

f2k — 1 2k + 1 

3. Outside the two sheets of the same cone lie hyperbolic lifts. For ^ ^ ( — ^ — ? — ^ — 

k an even integer, the points (p, a, 0) G -D x M, represent regular hyperbolic lifts, while 
for k odd integer they represent inverse hyperbolic lifts. We denoted by T-L2j, respectively 
^2j+i) J ^ ^ the subsets of regular, respectively inverse hyperbolic lifts. 

Denoting by d{£i, 'Hj) the boundary between two nearby regions we note that: 



4. The cones P. 
conjugated to 



1,1 

2k 
1 1 

1 



d{£2k-i,'H2k) \ {Tk}, A; G Z, contain lifts of matrices SL{2, 



i.e. matrices with b — c>0. The boundary P, 



1,-1 

2k 



d{£2k, T^2k) \ 



{Tfc}, A; G Z, consists in lifts G, whose projections are S'L(2, R)-conjugated to 



i.e. matrices with 6 — c < 0. 



1 




5. The cones V. 



-1,1 

2k~l 



d{£2k-iiH'2k_i) \ {T{2k-i)/2}i k E Z, contain lifts of matri- 



ces S'L(2, ]R)-conjugated to 



-1 1 
-1 



i.e. matrices with b — c > 0, while V. 



-1,-1 

2k-l 
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d{S2k,'H2k+i) \ {^(2fc-i)/2}5 k E Z, contains lifts of matrices S'L(2, M)-conjugated to 
-1 -1 



-1 



i.e. matrices with b — c < 0. 



In the study of twist numbers we are interested in the range, [m{G), M{G)], of the 
displacement map \1/ of a lift G. In order to deduce the position of such intervals 

for different regions in 5*1/(2, M), we show how are related the intervals [m{G), M{G)], 

[m{H-^GH), M{H-^GH)], associated to 5L(2^)-conjugated lifts. 

Proposition B.l The minimum, respectively the maximum displacements, of two SL{2, M 
conjugated homeomorphisms , belong to the same demi-unit interval, i.e. an interval of 
the form {k/2, [k + l)/2), for some k G Z.- 

Proof: 

We give the proof only for the minimum displacement of two SL{2, M)-conjugated 
homeomorphisms. More precisely we show that: 

T) Tl T) T) 

maxi- eZ/2\- <m{H-^GH)} = max{- G Z/2 | -< m(G)} 
2 2 2 2 

(44) 

min{- G Z/2 I - > m{H-^GH)} = min{- G Z/2 | - > m{G)} 



In order to get the first relation in (jH]) we prove that for any G G SL{2, M) we have: 

Tl Tl Tl Tl 

max{^ G Z/2 I ^ < x' - x} = max{^ G Z/2 | ^ < G{x') - G{x)} (45) 

Let x,x' G R be such that n/2 < x' — x < {n + l)/2, and distinguish two cases: i) 
x' = X + n/2, and ii) x' = x + 5, n/2 < 5 < (n + l)/2, ri G Z. In the former case we have 
G{x + n/2) = G{x) + n/2, i.e. G(x') — G(x) = n/2, while in the latter one, taking into 
account that G is increasing, we get: 

^(x + n/2) < G{x + 6)< G{x + (n + l)/2), i.e. n/2 < G(x') - G{x) < (n + l)/2 

In relation (HSl) for the liftiJ"^: 



maxi- G Z/2 I - < x' - x| = maxi- G Z/2 I - < H'^(x') - E~^(x)\, 
^2 ' ' 2 ~ ^ ^2 ' ' 2 ~ ^ ' ^ 

we take y = H^^{x), and x' = G{x) = GH{y), and get: 

max{^ G Z/2|^ < G{H{y)) - H{y)} = max{- G Z/2 | - < H-'GH{y) - y}, 

V y G M. Taking the minimum for ?/ G M we get the first relation in (jH]). 
Similarly one proves the second one. 

Because each region £k, 'H2k,'H2k_i, A; G Z is foliated by conjugacy classes, we have: 
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Proposition B.2 The displacement intervals for maps G in different regions of the 
S'L(2,M) are located as follows: 

aJIfGeSk, [m(G'),M(G)] C (A;/2, (A; + l)/2); 

bjIfGe U2k, then k G (m(G),M(G)); 

2k -\- 1 

c)IfGe H'^,^„ then G (m(G), M(G)); 

Proof, a) A conjugacy class in E^. is represented by tlie sliift T^^ witli uj G (A;/2, (fc + 
l)/2). The displacement map of is the constant function = w, and thus m{T^ = 
M{T^) G (fc/2, {k + l)/2). From Proposition ED we get a). 

b) A conjugacy class in 'H2k is represented by a lift, G, of the normal form, (HTj) . 
with m(G') = A; - 5, M(G) = k + 5, for some < 5 < 1/4. 

c) A conjugacy class in T-i^k-n k E Z is represented by a lift G of the normal H\ and 
m{G) = {2k - l)/2 - 5, M(G) = {2k - l)/2 + 5, for some < 5 < 1/4. 

Taking into account the dynamical properties of circle maps [17] , the properties of the 
translation numbers stated in|Xl and that the translation number of a shift is u we get 

that, if G G SL{2, M) is a lift of a circle homeomorphism g^, then the translation number 
of G is (see also [23], [ID]): 

1. an integer number, if G is a regular hyperbolic or a positive parabolic lift; 

2. a demi-integer, i.e. a number in (2Z+l)/2, if G is an inverse hyperbolic or a negative 
parabolic lift; 

3. a real number different from integers and demi-integers, if G is an elliptic lift. More 
precisely t(G) = 9 + k, k E Z, where 9 G (—1/2,0) or (0, 1/2), depending on the 
sign of b in n(G). 

In Fig{T3] a) we illustrate the values of translation number function, in each region in 
the space SL{2,R). 
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